心理生理交互分析(PsychoPhysiological Interactions,PPI)是一种用于分析任务态功能性磁共振成像(fMRI)数据的方法。它主要用于研究任务条件如何影响不同脑区之间的功能连接强度。
PPI分析的基本成分:
通过分析这三个项之间的关系,PPI可以揭示在不同任务条件下,脑区之间的功能连接是否发生变化。
PPI分析在SPM(Statistical Parametric Mapping)软件中有详细的实现步骤,可以帮助研究者进行任务态fMRI数据的功能连接分析。
PPI分析的步骤包括:
其中,想要生成PPI变量文件,需要我们使用一阶处理后生成的SPM,mat文件,二阶组水平上的SPM.mat文件由于缺少sess字段,不能够直接用于从零开始的PPI分析。所以之后的步骤我们以fz条件下0-6秒的单个被试的一阶SPM.mat为基础,进行一次PPI分析的流程(如图1)。
在我们完成一阶处理后,得到了个体水平上、不同condition的对比,这些对比写在了一阶生成的SPM.mat中,而选取VOI种子点就需要在这个基础上进行。这个过程我们可以想象是把原料(咖啡豆磨粉)放入加热壶得到加工产物(咖啡液体)的过程,咖啡豆磨粉是一阶SPM.mat、半成品咖啡液体是得到的VOI_*.mat。
首先在spm界面我们选择result进入结果查看,导入被试1_fz0-6_1st的文件,这里我们多项选择inhibition、distraction、distance仨condition,然后选择“done”(如图2)。
在之后的结果校正选项中,我们进行如下选择:conj'n、应用mask上我们选择none、多重比较校正选择FWE(p值设定为0.05)、阈值体素(Extend Threshold Voxel)设定为10(见图3)。进行好以上设定后得到统计参数图。
得到参数统计图后,我们仅仅是得到坐标,但是为了确定图中信号增强确切来说是哪个脑区,我们需要使用AAL3内置的MNI模板对坐标来打标签(如图4)。
打完标签后很明显,我们可以从参数统计图看到左右颞叶处存在两大块bold信号。我们假设要选取右边那个Temp_Sup_R作为VOI,那么首先我们可以点击一下感兴趣区域在列表左边的坐标系数,之后点击eigenvariate来尝试生成VOI,之后按照程序给出的参数进行填写(如文件名、VOI的形状、半径等等)。和之前在预处理中的生成文件一样,这里spm会在你的文件命名之前自动加上一个VOI_*前缀。
设置完毕后,我们将在原一阶SPM.mat路径处得到一个VOI_*.mat文件,这是Temp_Sup_R处Bold信号的时间序列。至此,我们拿到了“原料”之一–生理项素材。生理项代表的是特定脑区的BOLD信号变化。
进行PPI分析需要我们准备好生理项、心理项和交互项,其中在之前的步骤中我们已经准备好了生理项的文件载体,即右侧颞上回的种子点文件VOI_Temp_Sup_R.mat,而在这一部分的操作中,我们将生成心理项和交互项的文件载体–PPI变量文件。
通俗来理解的话,我们可以想象这一个步骤是将牛奶和咖啡倒入机器制作卡布奇诺的过程,而这里的牛奶和咖啡可以对应我们一阶SPM.mat和刚刚生成的VOI_*.mat文件。
在这一步骤中,我们可以直接从spm的主界面进入有关操作。我们可以直接在主界面点击PPIs选项,之后会直接蹦出对话框提示我们选择采用需要的一阶SPM.mat,这时候我们再选中之前说的fz条件下0-6秒的单个被试的一阶SPM.mat(见图6)。
选中后,程序会询问我们要进行心理生理交互(psycho-physiological interaction)分析还是生理生理交互分析,在选择psycho-physiological interaction之后,程序紧接着会弹出另一对话框提示我们选择对应的VOI_*.mat文件。刚刚我们说过,我们这次假定感兴趣的种子点是右侧颞上回,因而我们这里要选择的种子点文件便是在上一步生成的VOI_Temp_Sup_R.mat(如图7)。
这之后,程序会需要我们选择纳入的conditon和其在多元线性回归中的权重,这取决于我们之后的数据分析需求和实验设计以及研究目的,此处我们仅仅针对操作的流程,假定我们想要比较3strategies > nature这个condition,那么其中我们可以这样设置权重:
将三种策略inhibition、distraction、distance的权重分别设置为1
nature的权重设置为-3
blank我们选择不纳入回归方程。
注意,该参数之后再具体工作中的设置需要根据教授的要求具体情况具体修改。最后程序会让我们设置PPI变量文件的名称,这里我们假如输入“TemSupR_3vsNature”,程序也是会自动加上一个PPI_*前缀,因而我们最后得到的文件是PPI_TemSupR_3vsNature.mat。设置好后自动返回给我们一张结果图片。这里按照笔者标示,表(1)是VOI生理项Bold信号的时间序列;表(2)是心理项的信号时间序列;表(3)是心理生理交互项的信号事件序列(见图8)。之后,再VOI_Temp_Sup_R.mat的所在路径我们得到了心理项和交互项以及生理项的文件载体PPI_TemSupR_3vsNature.mat。
在之前的步骤中,我们知道PPI分析是将心理项、生理项、和交互项对目标脑区Blod信号时间序列来一个多元线性回归。而现在我们已经得到了VOI文件和PPI文件(即VOI_Temp_Sup_R.mat和PPI_TempSupR_3vsNature.mat)是不是就可以在回个归了呢。并不,我们之前总是比如说PPI_TempSupR_3vsNature.mat是心理项、生理项和交互项的文件载体云云。也就是说,我们想要对不同的项目来回个归还需要把其对应的值从文件中提取出来。这部分就是将心理项、生理项以及交互项的value从文件中取出来。
在这个步骤,我们可以首先创建一个n*3的Excel表格,其中n为你实际分析的图片数量或者说frame数量。而这三列依次可以是PPI交互项、VOI的Blod信号(生理项)、心理任务向量(心理项):
| 序号 | PPI交互项 | VOI生理项 | 心理项 |
|---|---|---|---|
| 1 | 数据 | 数据 | 数据 |
| 2 | 数据 | 数据 | 数据 |
| …. | …. | …. | …. |
| n | 数据 | 数据 | 数据 |
做成之后,我们接着就需要把实际的各项数据从载体文件中提取出来,复制粘贴进我们的整理表格了。这里假如我们首先去找生理项的数据。我们在MATLAB中打开对应的文件工作路径,加载VOI对应的mat文件:load('PPI_TempSupR_3vsNature.mat'),之后我们可以在工作区中看到该文件中有一个名为PPI的变量,我们查阅这个变量其中又包含Y、P和ppi这仨变量表格,它仨分别对应了生理项、任务项和交互项的值,把这仨变量表格里的数值复制粘贴到我们的Excel表格里备用就可以了。
当然其实我们也可以使用简单的代码命令来做:
% 加载PPI.mat文件
load('PPI_TempSupR_3vsNature.mat');
% 提取心理任务向量
Psych_variable = PPI.P;
% 提取PPI交互项
PPI_interaction = PPI.ppi;
% 提取从VOI提取的BOLD信号
V2_BOLD = PPI.Y;
在前文中我们提到过,PPI分析可以被看作多元线性模型在fMRI数据中的应用,那么既然是回归,我们便需要再次定义一般线性模型。而在spm的model specification模块中,便对应了我们需要把心理项、生理项和交互项放进回归方程里。这里的操作如果从0开始做的话整个流程和我们之前在培训时的数据一阶处理内容很像,这里其实我们可以直接将当时一阶的对应条件batch拿来改一改再用。
我们可以加载出之前fz0-6用的一阶batch,其他的条件和参数没有特殊要求可以不改,这里PPI分析只是在model specification中添加心理项、生理项、交互项这三者的回归(如图10),每个变量的名字name可以根据需要更改,对应的value就是刚刚我们在PPI_TempSupR_3vsNature.mat里往Excel表格里复制粘贴的那些值,这里对应填写上即可。
Regressor中添加PPI变量
改装完一阶batch的model specification后,我们还要去将其的contrast manager进行一定的小修改,具体而言就是添加一个新的t对比(t-contrast)。但是如果我们想要去查看脑区之间的连接,我们则需要添加交互项的对比项目,其重要意义如下:
PPI.ppi:显示了在不同任务条件下,VOI与其他脑区之间的功能连接如何变化。
PPI.P:显示了任务条件对脑区活动的直接影响,但不涉及功能连接。
PPI.Y:显示了VOI的时间序列,但在交互项中不涉及任务条件的影响。
不过这里有个疑问需要解决,我本来的VOI就是在特定任务条件下选取的,怎么就这里说不涉及任务条件的影响了呢?这是因为尽管原始的时间序列虽然的确反映了某区域在整个试验期间的活动,但是在生成PPI文件的时候这个时间序列会经过去卷积和再卷积处理,以去除与任务时间的直接关系。因而说PPI的这个生理Bold和任务没关系。
那么说到这里,我们继续进行最后一步的说明,即我们将PPI交互项加入contrast manager中。其实就是创建一个新的t-contrast,name就是写交互项比如“inter—PPI”、weight vector我因为之前定义回归模型的顺序是心理项、生理项和交互项,那么这里我的权重就是填[0,0,1](如图11)。
运行batch之后,我们得到了一个类似于一阶结果的SPM.mat同样我们使用spm主界面的result选项,查看交互项的对比,就可以观察哪些脑区的信号增强和我们VOI有所关联了。不够根据之前笔者观察,普通一阶结果和PPI结果在一些对比中会有信号显著性差异,经过询问ai得到以下,仅供参考:
在普通一阶分析和PPI分析中,信号显著性的差异可能源于以下几个原因:
- 分析目标不同
普通一阶分析:主要关注单个脑区在不同任务条件下的激活情况。它评估的是每个脑区在不同条件下的活动水平。
PPI分析:关注的是脑区之间的功能连接变化,特别是感兴趣区域(VOI)与其他脑区之间的交互作用。它评估的是任务条件如何调节这些连接 。
- 模型复杂度
普通一阶分析:模型相对简单,主要包含任务条件的回归因子。
PPI分析:模型更复杂,包含了任务条件、VOI的时间序列以及它们的交互项。这种复杂性可能揭示出普通一阶分析中未能检测到的信号。
- 信号来源
普通一阶分析:仅考虑单个脑区的活动。
PPI分析:考虑了脑区之间的交互作用,可能揭示出由于这些交互作用而产生的显著信号。
- 统计功效
普通一阶分析:可能由于统计功效不足,未能检测到某些脑区的显著信号。
PPI分析:通过考虑交互作用,可能提高了检测某些脑区信号的统计功效。